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Abstract The 4-dimensionally covariant approach to multiconstituent Newtonian 
fluid dynamics presented in the preceding article of this series is developed by con- 
struction of the relevant 4- dimensional stress energy tensor whose conservation in the 
non-dissipative variational case is shown to be interpretable as a Noether identity of 
the Milne spacetime structure. The formalism is illustrated by the application to ho- 
mogeneously expanding cosmological models, for which appropriately generalised local 
Bernouilli constants are constructed. Another application is to the lordanski type gen- 
eralisation of the Joukowski formula for the Magnus force on a vortex. Finally, at 
a global level, a new (formally simpler but more generally applicable) version of the 
"virial theorem" is obtained for multiconsituent - neutron or other - fluid star models 
as a special case within an extensive category of formulae whereby the time evolution 
of variously weighted mass moment integrals is determined by corresponding space in- 
tegrals of stress tensor components, with the implication that all such stress integrals 
must vanish for any stationary equilibrium configuration. 



1. Introduction 

This article is the second of a series providing a systematic treatment of the essential 
dynamical properties of multifluid models (as needed for the description of a neutron 
superfluid moving with respect to a normal background) using a non-relativistic but 
4-dimensionally covariant formalism in which the preferred Newtonian time gradient 
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features as a null eigenvector of the degenerate contravariant space metric tensor 7'*''. 
In the preceding aricle [1], which will be referred to simply as (1), it was shown how, in 
the absence of dissipation, the equations of motion are obtainable, in the conservative 
limit, from a variation principle given, in a fixed background gravitational field 0, by a 
Lagrangian A that is specified as a function of a set of conserved current 4- vectors n^. 
Subject to the condition that the variation of these currents be restricted to have an 
appropriate world line dependent form, the dynamical equations were derived as the 
condition of vanishing of a corresponding set of 4-force densities that were specified 
by the equation (1-159) in terms of the dynamical conjugates of the currents, namely 
the set of 4-covectors tt^ representing the effective energy-momentum per particle of 
the species labelled by the suffix X. 

The present article (II) describes the generalisation needed to allow for active self 
gravitational effects, which can be dealt with in the conservative case simply by adding 
an appropriate gravitational field contribution A^^j to the action density. It will be 
shown how the associated stress-momentum- energy density 4-tensor will be subject 
to various (generalised Bernouilli and virial type) conservation laws ensuing from the 
invariance of the theory with respect not just to Galilean transformations but also to 
their accelerated, Milne type, generalisations, (which are important for cosmological 
applications) . A formal derivation of the relevant Noether identities is provided in an 
appendix. A following article (III) will describe the generafisation needed to allow for 
non dissipative effects such as viscosity, resistivity, and drag on superfiuid vortex lines. 

2. Stress momentum energy tensor 

Experience with the technically simpler relativistic case [2] (and previous work on the 
Newtonian case [3] ) suggests the utility at this stage of introducing the stress momentum 
energy density tensor 

n = E<7r^ + *5^ (1) 

in which the covectors vr^ arc the 4-momenta defined by (1-152) and ^ is the generalised 
pressure function defined by (1-148). The utility of this (chemically covariant) tensor 
derives from the property that its divergence can be seen to be given just by the sum 

= (2) 

of the relevant non gravitational force densities - which could be expected to cancel 
each other out, even in a non-conservative model, so long as it is isolated from external 
interactions - plus the relevant gravitational force density which can always be expected 
to be present, by the simple formula 

V^n = /.-pV,0, (3) 

a relation that is shown in the appendix to be interpretable as a Noether identity. 

It is to be remarked that the 3-momentum density vector given by (1-146) will be 
obtainable directly from (1) as the mixed space and time projection specified by 

il'' = tpTP.Y" ■ (4) 
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The gauge independent character of this 3-momentum density il'*, as manifested by 
its original specification (1-145) is to be contrasted with the highly gauge dependent 
character of the corresponding energy density current, namely 

^ -y/^e'^ , (5) 

and in particular of its time component, the ordinary energy density, which will be 
given by 

C/ = C/% = -Enx7rX-*- (6) 
The latter will evidently be decomposable in the form 

U = U^., + U^^, , (7) 

in which the potential energy contribution is just the opposite of the corresponding 
Lagrangian contribution (1-94), i.e. 

t^pot = P0 = -Apot , (8) 
while the material contribution will be given by 

C/.at = -En,//^e'^-*. (9) 

The complete stress energy momentum density tensor will evidently have a cor- 
responding decomposition as the sum of a purely material part and a gravitational 
contribution in the form 

n = 3:/, - <i>P'u , (10) 

with 

T/. = En>- + *5^ (11) 

in which ^ is the same as before, i.e. there is no need to distinguish between ^ and 
^mat ' because its gravitational contribution cancels out so that it is directly expressible 
as 

* = A.at-EnX., (12) 
and its general variation (1-151) simplifies as 

S^^-EKSfxl. (13) 

Actually we can take the decomposition one step further by expressing the material 
contribution as the sum of a kinetic contribution and a purely internal (and therefore 
gauge invariant as well as chemically covariant) internal part in the form 

rp H _rp n .rp /-|^\ 

where the kinetic contribution is given simply by 

tJ.-T.KpI (15) 
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and the internal contribution is the Milne gauge invariant pressure tensor, P'^ say, 
which is given by 

Tf, = n = E ni^Xl + "^K , (16) 

with ^ given in terms of quantities obtained just from the internal contribution to the 
Lagrangian by the formula (1-148). The corresponding decomposition of the material 
energy contribution (9) will have the form 

U^..-U,,^ + U,..- (17) 
in which it can be seen, using (1-48), that we shall have 

t^Hn = = 7>xV. , - E Y.Kxl - ■ (18) 

In terms of the ether rest frame chemical potentials defined by 

X" = -X".e'' (19) 

it can be seen that we shall have 

= E n,x" - * , sa^^ = E + E ^xl ■ (20) 

It follows from the identity (1-144) that the contraction of the pressure tensor with 
the degenerate space metric ^'^'^ gives a result that is symmetric, i.e. 

7^[^P^j, = . (21) 

It can be seen that this entails a corresponding symmetry property for the space pro- 
jected part of the complete stress momentum energy density, namely 

7^M7.].n = 0- (22) 

Putting all the pieces together again, we see that this complete stress momentum 
energy density tensor will be expressible as 

n = n + E K(PI - m^ct>U) . (23) 



3. Action and stress energy for the gravitational field 

Up to this point we have been treating the gravitational field just as a given back- 
ground, but we can promote it to the status of an active dynamical field by taking to 
be an extra independent variable in the Lagrangian and adding in an extra gravitational 
field term so as to obtain a total action density given by 

A.„, = A + A,, , (24) 
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in which the gravitational field term is given by 

A^^, = (87rG)-VV^0, (25) 
where the graviational field is defined according to the potential as 

= -r''V.(t>- (26) 

Since the only other contribution involving </) is the potential energy term Ap^^ given 
by (1-136), it can immediately be seen that the requirement of invariance of the total 
action with respect to localised perturbations of does indeed give the usual Poisson 
source equation, which, in the covariant formulation we are using here, will be given by 

7'^'V^V,0 = 47rGp. (27) 

This can be alternatively presented in a manner reminiscent of the relativistic Einstein 
equation, using the notation of the Newton Cartan formalism described above, as 

R^^^AnGpt^t,, (28) 

where R^i, is the Ricci type trace of the Newton Cartan curvature, as given by (1-31). 

For the reason discussed in the appendix, there will be an associated gravitational 
stress energy momentum density tensor given by the formula 

T = -(47rG)"i(5^V.0 - \6^^9''V,<P) (29) 

(which was misstyped in the previous version [3] where the initial minus sign was omit- 
ted). 

The purely gravitational contribution (29) can be added to the contribution (23) 
associated with the material source distribution to give a grand total 

T "^^T^^ + T H, (30) 

that satisfies the identity 

v.3:„r. = /. , (31) 

in which all that remains on the right is the sum (2) of the non gravitational force 
density contributions, if any, that are defined by the variational formula (1-159). 
In the strictly conservative case for which the variational field equations 

/^ = 0, (32) 

are satisfied, and much more generally in any system that is effectively isolated so that 
the separate contributions in (32) cancel out in the sum (2), leaving a vanishing total 
force density, /^^ = 0, we shall simply be left with a total energy momentum conservation 
law of the form 

^.TJ. = , (33) 
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in which the conserved tensor field (30) can be redecomposed in the form 

r't = T'^ +T^,, (34) 

toti^ mati^ grtl^ ' V / 

where the total gravitational contribution is given by 

T = T ''^ + T . (35) 

It is to be noted that there will be a corresponding total energy density 

U... = -t.TJ.e'' + U^, (36) 

in which the contribution (29) provides a gravitational field energy density that can be 
evaluated, using (5), as 

t^g. = -t.Tjy - -K. = -\p^ - {STrGy'wMgn ■ (37) 

It can thereby be seen from (7) that the total energy density will be expressible in the 
form 

t^tot = U^.. + (38) 
in which the total gravitational energy contribution 

U^. = C^po. + C^g. (39) 

will be given by 

= -^M^^^r.e^ = \p<t> - {SnGr'V,{cf>gn . (40) 

It is to be observed that the final divergence term in (37) and (40) will cancel itself 
out when integrated over the volume surrounding a confined source provided that (as 
will always be possible in such a case) the gauge is chosen in such a way that the 
large distance limit of is zero. The remaining term in (37) will then half cancel the 
corresponding potential energy contribution given by (8), leaving a remaining term in 
(40) that is only half of U^^^. The effect of this semi-cancellation is to provide a total 
in which the gravitational binding between each pair of particles is (as it should be) 
counted only once, not twice. 

Another point to be noted is that there will be no corresponding gravitational field 
contribution to the space-like 3-momentum density (4), i.e. we shall have 

tpTJ^Y' = . (41) 
It can thus be seen from (5) that (4) will be replaceable by the equivalent specification 

= tpT,:.r' ■ (42) 



4. Space time symmetries and homogeneous cosmological models. 
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Let us now consider cases in which all or part of the system under study is invariant 
with respect to the action of one or several symmetry generators, k^, a = 0, 1, .... (For 
example in an application to a neutron star model, kQ might be a stationarity generator, 
i.e. the generator of a one parameter time translation group (in which case it would 
specify a natural choice of ether vector by the identification — k^) and /cf might 
be an axisymmetry generator, i.e. the generator of a one parameter group of rotations 
about some preferred symmetry axis. This means that all physically well defined fields 
must be invariant under the action of the corresponding Lie differentiation operators 
which we shall designate by the short hand notation £a = ka.£, or simply £ = k£, 
when we are only considering a single generator so that no ambiguity arises. It is to be 
emphasised however that physical invariance does not necessarily require vanishing of 
the Lie derivative of a field q that is gauge dependent, but only that its Lie derivative £g 
should be cancellable by some infinitesimal gauge transformation of the kind denoted 
by dq in the discussion at the end of Section 6 of [1]. More specifically, all that is 
required is that for each of the relevant Lie differentiation operators £a there should be 
a corresponding infinitesimal gauge transformation operator da such that any relevant 
field satisfies the condition 

£ag + dag = . (43) 

The occurrence of symmetries of this more general kind, involving non-vanishing in- 
finitesimal gauge transformations, is exemplified by the noteworthy case of the Milne 
type [4] homogeneous cosmological models that are described below. 

A minimal requirement for any symmetry generator k^ is that its action should 
preserve the Milne structure of Newtonian space time. For the basic Coriolis structure 
fields and 7^*^ the question of whether there is a corresponding gauge transformation 
operator da does not arise, since they are gauge independent, so all that is required is 
the vanishing of their Lie derivatives, conditions that reduce by (1-17) to the form 

uv^k:: = , y^^^pK^ = , (44) 

from which the space projected and complete space time divergence conditions 

7^VA' = 0, V.A;: = 0, (45) 

are obtained as corollaries. In a particular Aristotelian frame with respect to an ordinary 
orthonormal system of coordinates the conditions (44) will be expressible as 

\k^ = 0, Vik°=0, V(V^ = 0. (46) 

In the context of general relativity the well known condition for k^ to be a symmetry 
generator of the spacetime structure is just the well known Killing equation V^^kg^^ = 0, 
which guarantees the invariance of the spacetime metric gf^j^ and hence also of the 
derived Riemannian connection. However in the Newtonian case, since the connection 
is not simply derivable from the structure fields and 7^*^, the corresponding first 
order differential conditions (44) will not by themselves be sufficient to quahfy k^ as 
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a space time symmetry generator. In order for to qualify as a Newtonian space 
time symmetry generator its action must also preserve the fiat connection T^'^p- Since 
the latter is subject to the gauge dependence condition (1-75) its Lie derivative will 
not have to vanish absolutely, but only modulo the action of some infinitesimal gauge 
transformation, which will be characterised by some (infinitesimal) boost potential /3a, 
with a derived time (but not space) dependent boost vector and a corresponding 
acceleration vector a'^, that are given, as in (1-15) by 

= t'^'^ V./3a , t'^'^VA^'^O. (47) 

and 

< = e'^^^K = T^'^V.aa , «a = V,/5a . (48) 

Modulo a choice of rest frame at some arbitrary reference event, the specification, 
according to (1-16), of an ether frame field is equivalent to the specification of a 
corresponding connection F^'^^, namely the one with respect to which the vector field 
e'* is mapped onto itself by parallel transport. The preservation of this connection by 
the action generated by will therefore be ensured simply by the requirement that it 
should preserve the ether vector. Substituting this vector in place of q in the general 
purpose preservation condition (43), one sees from the transformation rule (1-74) and 
the parallel transport property (1-16) that the requirement of preservation of the ether 
vector is equivalent to the condition that the relevant boost transformation should 
be given simply by 

K = e^^.k^ . (49) 

It is evident that this will conveniently vanish, so that there will be no need to bother 
about allowance for the gauge adjustment, in cases for which the symmetry generator 
k^ is time independent, as will be the case for the symmetries that are most relevant 
in applications to rotating star models, namely stationarity and axisymmetry. On the 
other hand the gauge adjustment will have an indispensible role in the kind of time 
dependent translation symmetries that are relevant in homogeneous self gravitating 
(cosmological type) configurations such as will be described immediately below. In all 
cases, it can be seen from (44) that the formula (49) can be used to express the gradient 
of the symmetry generator as the sum of orthogonally projected and mixed components 
in the form 

vA^ = 7;7;vA"+i^6:. (50) 

In terms of the boost acceleration obtained from (49) using (48), the invariance 
condition for the connection F^^^^ will be given according to the general principle (43) 
by the formula 

V,V,k^ = t,Ua^, (51) 

whose derivation is based on the use of the simple fiat space special case of the Yano 
formula [5] for the Lie derivative of a symmetric connection. 

The general version of the Yano formula, including allowance for curvature, is needed 
for the application to the Newton-Cartan connection (1-25) whose Lie derivative will 
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have the form given by 

£,u/,^D,D,k<' + R,^\k'' . (52) 

In the Newtonian case (unlike the relativistic case) the conditions for invariance of the 
relevant (Milne) space time structure, namely the conditions (44) and (51), are not 
sufficient to ensure invariance of the gravitational field, as embodied in the independent 
gauge invariant connection a;^''^,, or equivalently in the gauge dependent field g^. In 
view of its gauge independence, the invariance under the action of of the Newton 
Cartan connection requires simply that its Lie derivative, as given above, should vanish, 

A<p = 0. (53) 

Subject to the Milne structure invariance conditions (44) and (51), it can be seen from 
the formula (52) (using the expression (1-29) for the curvature) that the supplementary 
condition (53) for invariance of the gravitational field reduces to the condition obtained 
by application of the general requirement (43) to the gauge dependent gravitational 
field vector g^. It evidently follows from (1-23) that this invariance requirement will 
take the form 

£-.9' = <. (54) 

in which the Lie derivative of the field will of course be given by the well known com- 
mutator formula 

£,g^ = k:V,g^-g^V,k^. (55) 

It similarly follows from (1-24) that at the more highly gauge dependent level of the 
gravitational potential 0, the corresponding field invariance condition will be expressible 
simply as 

k^^VA^-oc., (56) 

where eta is the acceleration potential given by (48). 

The prototype illustration of symmetries that are non-trivial, in the sense of re- 
quiring a non zero gauge transformation according to (49), is provided by the case of 
homogeneous isotropic cosmological type models, with expansion described in terms of 
a comoving radial scale factor a say, and 3- velocity given in terms of Cartesian space 
coordinates X'^ by v'^ — HX^ where the (time dependent but spacially uniform) Hubble 
parameter is given by H = &/a (using a dot for ordinary time derivation). Since the 
pressure is postulated to be uniform it has no effect on the motion, so the fluid particles 
will be effectively in free fall with respect to the background gravitational field, which 
will be given by = |7rGp7ijX*X-'. 

One of the things that delayed the discovery of these configurations is their lack of 
stationarity, i.e. their essentially time dependent nature. Their matter distributions, 
and in particular the density p are not invariant with respect to the action of the 
ordinary time translation symmetry generator kg that is definable in terms of the fiat 
Aristotelian coordinates {t, X*} used above by 

A;o° = l, A;o' = 0, (57) 
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which means that it is simply identifiable with the relevant Aristotelian ether vector, 

= . (58) 

The misguided presumption that the cosmological model would have to be stationary 
was one of the psychological obstacles that, prior to its demolition by Friedmann's theo- 
retical insight (and Hubble's observations) had prevented everyone (including Einstein) 
from seriously attempting to carry out the non trivial but (compared with the devel- 
opment of general relativity) not so difficult technical step that was finally taken by 
Milne. 

The stationarity generator (58) illustrates a distinction that does not arise in General 
relativity theory, where the property of being a Killing vector, i.e. a solution of the 
equation V^^k^^ — 0, ensures the invariance not only of the background metric g^y but 
also, automatically, of the gravitational field as embodied by the associated connection. 
In a Newtonian context it is necessary to make a distinction between what may be 
termed weak Killing vectors, or Killing-Milne vectors^ meaning those that just satisfy 
the conditions (44) and (51) for invariance of the background spacetime structure, 
and what may be termed strong Killing vectors or Killing- Cartan vectors, meaning 
those that not only satisfy the conditions (44) and (51) for preservation of the Milne 
structure but also the supplementary condition (54) for invariance of the gravitational 
field and hence of the complete Newton Cartan structure. Regardless of the matter 
distribution, the stationarity generator kg given by the specification (57) will always be 
a KiUing-Milne vector, since it obviously satisfies (44) and also (with vanishing boost 
acceleration) (51). However in the particular application described above k^ is not a 
Killing vector in the strong sense because it does not satisfy (54). 

Although they are not stationary, the cosmological configurations described above 
are obviously isotropic in the sense of being invariant (in this case without any need 
for an associated boost transformation) under the action of the set of ordinary angular 
symmetry generators kf that are definable for values j = 1, 2, 3 of the index a as follows. 
For a rotation about an axis specified by a space unit vector z/j', with components given, 
with respect to the fiat Aristotelian coordinate system we have been using, simply by 

u^ = S^, (59) 

the corresponding angular symmetry generator will be given by the prescription 

kP^O, = (60) 

in terms [6] of the fiat space Killing- Yano 2-form, whose (Cartesian) components are 
given by 

Yij^eijkX\ (61) 

and which is characterisable, even with respect to non Cartesian coordinates, by the 
Killing- Yano equation 

V(,y,), = 0. (62) 
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Whatever the matter distribution may be, it is evident that, Ukc the stationarity 
generator A;,^, these axisymmetry generators k^^ for j = 1,2,3 will always satisfy (44) 
and (again with vanishing boost acceleration) also (51), so that they will always qualify 
as KiUing vectors in the weak sense. In the particular example of the cosmological 
application described above it is also evident that (in this case unlike the stationarity 
generator ko") the rotation generators will satisfy the further condition (54) that qualifies 
them for description as Killing- Cartan vectors. 

We now come to the most essential, though not so obvious, symmetry feature of 
the cosmological configurations described above, which is that as well as being isotropic 
they are also homogeneous (with respect to space but not time) in the sense of being 
characterised by a set of time dependent space translation symmetries with generators to 
which we shall attribute negative label values, a=-l,-2,-3 (since the positive label values 
a=l,2,3 have already been used up in the specification (60) of the rotation symmetry 
generators). With respect to the same ordinary Aristotelian system of orthonormal 
coordinates {t, X^} as above, these space translation symmetry generators will be given 
by the specification 

k^. = , kl^^ aSj , (63) 

for j — 1, 2, 3. What delayed the discovery of such symmetries (even after Friedmann 
had overcome the psychological barrier of the stationarity presumption) until Milne 
saw how to exploit the gauge invariance, is that (in order for the generators (63) to 
quahfy as Killing vectors in both the weak and strong sense described above) the effect 
of their action needs to be cancelled by the effect of corresponding (non-hnearly) time 
dependent (so non Galilean) transformations characterised by (spacially uniform) boost 
vectors b2j and boost potentials given by the formulae 

b:^^a5i, /3_j=d7j,X\ (64) 

In particular the velocity field can immediately be seen to satisfy the invariance condi- 
tion £_a f * = 61^, while the gravitational field will satisfy the corresponding invariance 
conditions 

£_j^^ = a4, = (65) 

(so that the translation generators /cj^ qualify as Killing vectors in the strong sense) 
where 

a_*j = aS-' , Q;_j = (^JjiX'^ ■ (66) 

By contraction with the stress momentum energy density tensor T'^, we can use the 
symmetry generators k^ to construct corresponding generalised momentum currents 

= nk: , (67) 

of which, for example Vq^ will be interpretable as the negative of an ordinary energy 
current, Vi" as a current of ordinary angular momentum (about the X'^ axis in the 
Aristotelian coordinate system used above) while V^i will be interpretable as a current 
transporting a kind of generalised linear momentum (in the direction of the X'^ axis). 
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The latter would reduce to a current of linear momentum of the ordinary kind in the 
non expanding (weak gravitational coupling) limit characterised by a time independent 
value of the factor a in the specification (63). 

Using the Killing vector conditions (44) in the decomposition (50), it can be seen to 
follow from the stress momentum energy tensor symmetry property (22), that we shall 
have 

nV.fc: = T%tX , (68) 

where is the boost vector given by (49). As a corollary of this useful lemma, it evi- 
dently follows from the tensor divergence formula (3) that the corresponding generalised 
momentum densities will have scalar divergences given by 

"^^K = K{u - PV.0) + il'^V./Ja . (69) 

On the right of this formula, the first term (involving the total external force if any) 
is of a familiar kind, representing the ordinary rate of energy loss (per unit volume) in 
the case of the stationarity generator /cg* and representing torque density (about the 
axis) in the case of the axisymmetry generator /cf , and in these cases the last term 
will be absent. This last term (involving the gradient of the boost potential /3a) will 
only be needed in cases for which the symmetry generator is time dependent. In the 
case of an isolated system, as characterised by a vanishing total (external) force density, 
= 0, and by a total mass current that satisfies the Newtonian mass conservation 
law (1-108), it can be seen that (69) will be expressible in the form 

V,{V: - (5^) = -p(aa + KVA) , (70) 

in which the right hand side will drop out altogether if the gravitational potential 
satisfies the condition (56) of invariance under the symmetry action generated by k^, 
a condition for which it is necessary that should be not just a Killing-Milne vector, 
as has been assumed in the derivation of (69), but more particularly a Killing-Cartan 
vector in the sense discussed above, meaning that its action preserves the gravitational 
field as well as the background spacetime structure. 

Independently of this latter (field invariance) restriction, but subject to the require- 
ment that the boost contribution should be absent, an isolated system will always be 
characterised by a conservation law for a corresponding, suitably adjusted, total gener- 
alised momentum vector. To start with, we define the (unadjusted) total momentum 
vector in the obvious way, simply replacing T'^ in (67) by the total stress energy tensor 
(30) which gives 

= TJX. = + T^X ■ (71) 

In the case of an isolated system, as characterised by a vanishing total (external) force 
density, = 0, and by a total mass current that satisfies the Newtonian mass 
conservation law (1-108), it can be seen that the analogue of (69) for this total will 
reduce to a generalised total momentum divergence law of the form 

V.iV^: - (3y) = -pa,, . (72) 
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Using the Poisson source equation (27) and the identity a^Vyg^ = Vu{aa,g^ + (pa^) the 
terms can finally be regrouped in a single divergence on the left in the form 

^•^KZ = , (73) 

expressing conservation of a suitably augmented total generalised momentum current 
that is given by the prescription 

'P.Z = - P^P" - i^^G)-\a^g>^ + <Pa:) . (74) 

The augmentation contributed by the second, third, and fourth terms will actually be 
needed only in cases for which is time dependent (as in the cosmological example 
characterised by (63), (64) and (66) as described above). The first term is all that 
remains, i.e. the unadjusted total V^^^^'^ will be conserved as it stands, in the more 
familiar cases in which k^^ is a generator of ordinary time translations (in which case 
the conserved vector will be proportional to energy flux) or rotations (for which it will 
be interpretable as angular momentum flux). It is to be emphasised that the validity of 
the generalised total momentum conservation law (73) follows just from the postulate 
that k^ is a Killing-Milne vector in the weak sense discussed above (meaning just that 
its action preserves the background spacetime structure), so that (unlike the results to 
be derived in the next section) it will be valid regardless of whether or not the relevant 
material field configurations have any corresponding symmetry property. 

5. Fluid symmetries and generalised Bernoulli theorems 

In contexts for which the local gravitational coupling can be considered to be suffi- 
ciently weak, it can be of interest to consider configurations of the generic kind governed 
by (69) in which a material medium, such as the multiconstituent fiuid dealt with here, 
is not necessarily subject to the symmetries of the background space time structure and 
the background gravitational field. 

The purpose of the present section is to show how much stronger conclusions can 
be drawn if the medium is itself invariant under the action of a Killing Cartan vector 
k^ of the kind characterised by the space time symmetry conditions (44), (45) and the 
gravitational symmetry condition (54) discussed above. 

In the case of a multiconstituent fiuid, it can be seen from the formula (1-77) for 
dTT^ that in terms of the boost potential, the corresponding symmetry conditions on 
the momentum covectors will be given, according to the general principle (43), by 

£a7r^ =m^V./3a. (75) 
Since, by its definition, the Lie derivative of the momentum covector will be given 

by 

£^nl = 2A;:V[,7r-^] + V^(A;:7r\) , (76) 

we see from the definition (1-161) of the vorticity that the symmetry condition on the 
momentum will be expressible in the form 

= V^B,- , (77) 
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in which is a generahsed momentum component that is defined by the formula 

B^ = k:nl-m^f3,, (78) 

and that can be identified as a generafisation of the historic Bernoulh scalar. 

It is immediately obvious from (77) that we obtain a strong generalisation of the 
Bernoulli theorem to effect that if the ffow of some particular constituent is irrotational, 
then the corresponding Bernoulli scalar will be uniform throughout, in the sense of being 
independent of both space and time, 

^^.-0 ^ V.i3- = 0. (79) 

In the case (of the kind originally considered by Bernoulli) of ordinary time transla- 
tion symmetry, as obtained by setting a=0 with the specification (57), the corresponding 
value, 

= i^le^ , (80) 

will be interpretable as the negative of an effective energy per particle. Similarly in the 
case of an ordinary rotation symmetry, as obtained, for example by setting a=l (for a 
rotation about the X'^ axis) according to the specification (60) the corresponding value, 

= TT-fcr , (81) 

will be interpretable as an effective angular momentum per particle. In the case of a 
boosted translation symmetry along for example the obtained by setting 

a = — 1 in the specification (63), we obtain a generalised Bernoulli scalar 

B^^ = k^^T^l - m^f3_, , (82) 

of a new kind, whose interpretation is not so elementary, due to the involvement of the 
boost potential 

In the general case, involving rotation or even transfer of matter between different 
constituents, it can be seen by combining (77) with the general force density formula 
(1-164) that the gradient of the relevant generalised Bernoulli scalar along the flow lines 
will be given by 

ny.B- = Kf- + KnlV^ . (83) 

We thus obtain a generahsed weak Bernoulli theorem to the effect that in the conser- 
vative case of vanishing force density, the Bernoulli scalar will be constant along each 
separate flow line, 

fl^Q, V^^Q ^ n^W^B^^O. (84) 

A useful alternative presentation of this generalised Bernoulli theorem can be ob- 
tained by introducing generalised momentum current densities that are defined for each 
separate constituent by the prescription 

Vr-^IK- (85) 
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The theorem (83) can thereby be reformulated as a divergence law of the form 



V 



(86) 



It can be seen that the generalised total (non gravitational) momentum density (67) 
is related to the sum of the separate constituent contributions by 



It was already made apparent by (71) in the preceeding section that the combination on 
the left of this equation will be conserved under rather general circumstances: all that 
is required is that the system as a whole should be effectively isolated from external 
influences (other than gravity) and that the generator should be a Killing- Cartan 
vector in the sense specified above. What has been shown in the present section is that if 
we make the more restrictive postulate that the material system itself is invariant, in the 
gauge adjusted sense characterised by (43), under the action generated by the Killing- 
Cartan vector, and if we also postulate that the separate force density contributions 
(and hence also the separate decay rates Vy.) all vanish individually, then (86) will 
be applicable so as to provide the much stronger conclusion that each of the distinct 
contributions in the sum on the right of (87) will be separately conserved. 

6. Generalised Joukowski theorem 

As an application of the generalised Bernouilli theorem discussed in the previous 
section, the Joukowski formula [7] is extended to the case of multiconstituent fluid. We 
shall derive the Magnus force arising on a perturbing vortex moving in an asymptoti- 
cally uniform medium characterised by vanishing currents n^^ = (asymptotic values 
will be denoted by an overhead bar). The gravitational potential is supposed to be un- 
affected by the perturbation. We shall follow closely a previous analysis in a relativistic 
context [8]. The fluid flow is supposed to be stationary thereby admitting the ether 
vector as a Cartan Killing vector, and longitudinally invariant along the uniform 
space like Cartan Killing vector l^ (vortex symmetry axis) satisfying IH^ = 0. 

The multiconstituent flow is described by the conservation law (33) V^T^^'j, = 0. 
Besides each of the currents are separetely conserved 'V^n^ — 0. It is further assumed 
that not only the total force density vanishes, but that each fluid component is isolated 

— n^va\^ — 0. It thus implies that the generalised vorticity 2-form will vanish 
whenever it is initially zero = 0. The framework in which the Joukowski formula 

is derived, is restricted to such irrotational flows (whether it is superfluid ot not). 

Consequently, from (79), the Bernouilli scalars associated with the corresponding 
Killing vectors are uniform Vi,Bl — 0. In particular, the following BernouiUi scalars 
are constants: 



(87) 



fix _ ^Xg. = ^^g^^ 



(89) 



8. 



;8 
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The multiconstituent fluid is exerting a force (per unit length) on the vortex given 

by 

^.-fTJ^'e^.dx^, (90) 
in which the "background" antisymmetric measure tensor is defined by 

The conservation of the total stress energy momentum tensor allows one to evaluate 
the integral along any closed circuit, which for convenience can be chosen sufficently 
far away from the vortex core for a linear expansion to be valid: 

T,:. + STJ^ + 0{S') . (92) 

The force is zero by definition in the unperturbed uniform medium hence it reduces to 

jr. = + 0(5") . (93) 

The linear deviation of the stress energy momentum tensor from the uniform back- 
ground value is expressible as 

STJ^ = E (^. Sn^ + n^^ Sn^ - n^^ Stt- S^) . (94) 

The force per unit length to lowest order is therefore given by SJ^i, — § ST^^^^ ^e^ndx'^ , 
i.e. 

5Tu ^T^i^.j Sn^ *e^i,dx^ + 2n^^ ^ Sn^ X]^.dx''^ . (95) 

It is then worthwhile to notice as a consequence of the Bernouilli theorem that the 
variation of the momentum is purely orthogonal to the vortex 

±1 57r- = Snl , (96) 

where ±^ is the operator of projection orthogonal to the vortex (i.e. orthogonal to 
the Killing vectors). The spacetime metric 61 is decomposible into the sum of the 
operators of projection parallel rj^ and orthogonal ±^ to the vortex. The contravariant 
antisymmetric measure tensor *e'^'' is introduced by 

= tpl^e"""" . (97) 

The covariant component l^r of the KiUing vector is well defined since this vector is 
spacelike. The projection operators are given by 

< = e'^U + riu (98) 

(99) 
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The force per unit length therefore becomes 



= E (j^u j K "e^^dx^' + 2n^" j ^tt^ 'e^^^dx^^ . (100) 
Using the following identity, 

(101) 

the force per unit length acting on the vortex simplifies to 

(5J^. = E(^.(^^x + ^"*£a.5C^), (102) 

in which the momentum circulation integral and the current outflux integral "D^ are 
defined by 

^ f K*e,^dx^ (103) 



= y TT^ dx" . (104) 

The current and vorticity conservation laws ensure that and respectively, do 
not depend on the integration closed path. Therefore the circuit can be chosen to lie 
at a very large distance from the vortex core so that the force is exactely given by the 
linear perturbation term. In cases for which there is no current creation in the vortex 
core, the current outflux integral will simply vanish V^. — and since the asymptotic 
values of these two integrals must also be equal to zero the force per unit length acting 
on the vortex is eventually given by 

^.=*£a.E<^''^''- (105) 

7. Virial moment theorems for isolated system 

The preceeding analysis has been essentially local, but of course whenever one is 
concerned with a system - such as a star model - that is confined within a compact 
region it is of particular interest to consider global quantities, particularly those that 
are subject to simple evolution equations or actually conserved. 

A specially noteworthy example is of course the total energy E^^^, which is definable 
as a function of time by a space integral that will be expressible in the form 



tot 



= / t^tot d'X , (106) 



with respect to standard Aristotelian coordinates {X",X*} of the kind described in 
Section 2 of [1], with time coordinate X" = t, and space coordinates X* of ordinary 
Cartesian type, so that the 3-dimensional Euclidean metric components are those 
of a 3 by 3 unit matrix, while the time covector has components — 1, ti — 0, and the 
corresponding ether velocity vector has components e° = 1, e* = 0. 
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It is evident from the local energy momentum conservation law (33) that the time 
evolution of this quantity will be expressible as the space integral of a divergence in the 
form 

^£^.o. = /v.TJ„d3X. (107) 

Using Green's theorem to convert this to an asymptotic surface integral at large dis- 
tance, it can be seen to follow that for an isolated system the energy will actually be 
conserved, i.e. we shall simply have 

^^.o.=0, (108) 

provided that (as is always possible in such a case) the gauge for the gravitational 
potential in the formula (29) for the external field contribution TJ^^ is chosen in 
accordance with the usual convention that it should tend to zero at large distances. 
(This condition automatically entails that will fall off as the inverse of the radial 
distance, a requirement which is sufficient to get rid of the boundary term that would 
otherwise be left over.) 

The foregoing total energy conservation law depends only on the time projected 
part of the complete local energy momentum conservation law (33), which can be 
decomposed in terms of separate energy density and 3-momentum density components 

u _rpo n' = T°' = y'^T % , (109) 

tot totO ' tot ' totfC ' ^ ' 

in the form 

Vof/.„. = V.T^:„, (110) 
(the part from which (107) is derived) and 

= -v,Tj; = -y^v.Tj, . (Ill) 

(in which, as remarked above there is no need to distinguish between 77* and 77 * since 

V ' *^ tot 

there is no purely gravitational 3-momentum contribution). The space projected part, 
namely the local 3-momentum conservation law (111) is particularly informative when 
used in conjunction with the conservation law (1-108) for the Newtonian mass current, 
whose time and space components 

p'^W (112) 

must satisfy 

V,p^-Vip\ (113) 
By combining this with (111) we obtain a noteworthy second order differential relation, 

Vo v„p = v.v,^J^ (114) 

that may appropriately be referred to as the local Newtonian virial equation, since, 
as will be shown below, it is what ultimately accounts for various older (specialised) 



18 



and newer (more general) variants of what is commonly referred to by the term "virial 
theorem" . 

The relevant global applications of the foregoing local conservation laws involve 
global mass - moment integrals of the form 

J pXd^X (115) 

and momentum - moment integrals of the form 

I wy^d^x (116) 

where X and are weighting factors that are constructed as fixed (i.e. time inde- 
pendent) given functions of the Cartesian space coordinates X*. In the same way as 
(108) was obtained by applying (110) to (106) using Green's theorem to get rid of a 
divergence, it can be seen that - for a locally confined system — application of (113) 
in (115) gives 

M = I WV^Xd^X. (117) 



dt 

In order for application of (111) in (116) to provide the analogous relation 

it is not enough to require that the material system be locally confined: in order for 
Green's theorem to be able to get rid of the relevant divergence contribution the weight- 
ing factor yi must satisfy an appropriate admissibility restriction. Since, according to 
(29), the components of the gravitational field contribution T*] will fall off as the fourth 
power of the radial distance from the material source, it can be seen that the criterion 
for admissibility in (118) is that the components y^ should grow more slowly than the 
square of the radial distance. Subject to the proviso that X should satisfy this ad- 
missibility condition when one makes the substitution y^ = ViX, it can be seen by 
combining (117) and (118)- or directly from the local virial relation (114) - that we 
shall obtain a generic global virial relation of the form 

jTJV,V,Xd'X. (119) 

The simplest example of the application of the foregoing formulae is of course to the 
case of the total mass M, which is obtained simply by taking the weighting factor X 
to be unity. Thus by setting 

X^l =^ M^M, ViX^O, (120) 

we see that (115) and (117) reduce respectively to the definition and conservation law 
for the total mass as given by 



M 



jpd'X, f^M^O. (121) 
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The next simplest possibility is the case of the dipolc moment, which is obtained by 
taking the weighting factor to be linearly dependent on the Cartesian space coordinates 
X'^. Thus in particular the dipole moment in the direction of let us say the axis is 
obtained just by taking the weighting factor X to be X^. Thus by setting 

X^X^ =^ M^V\ ViX^Yi: ViVjX^O, (122) 

we see that (115) and (117) reduce respectively to the definition and variation law for 
the dipole moment component as given by 

p3 ^ j pX' d^X , ^ J ' 

where the latter integral is interpretable as the ordinary linear 3-momentum in the 
direction of the X^ axis, whose conservation is now given by the relevant application of 
(118), or equivalently of (119) which can be seen to reduce simply to 

^P3^0. (124) 

This particular application will of course reduce to a triviality if we exploit the freedom 
to use the centre of mass frame characterised by the condition that the dipole moment 
simply vanishes. 

A less trivial application is to the case of angular momentum, whose components are 
obtained by taking the weighting factor to be given by the corresponding components 
of the Killing - Yano 2-form (61). Thus in particular the angular momentum Jg about 
the X^ axis will be obtained by identifying the weighting components yi with the Killing 
- Yano components Yi^ — Si^kX'^ ■ Thus by setting 

yi^Yi^ ^ j^j,, Viyj^Eij,, (125) 

we see that (116) and (118) reduce respectively to the definition and the conservation 
law for this angular momentum component, as given by 

J3 = e,j, I X' Wd'X , ^-^3 = 0, (126) 

in which the last step, namely the vanishing of the time derivative, results from the 
stress tensor symmetry property 

ri^] = 0, (127) 

which holds as a consequence of the corresponding property (22) of the material con- 
tribution, and of the manifest symmetry of the gravitational contribution (29). The 
familiar conclusion that the total angular momentum of an isolated Newtonian system 
will be conserved can alternatively be interpreted as an immediate global consequence, 
via Green's theorem, of the relevant particular application of the general local conser- 
vation law (74). 
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Having seen how mass and angular momentum conservation laws of the usual kind 
can immediately be recovered as particularly simple special cases for which the weight- 
ing factor is uniform or linearly dependent on the Cartesian space coordinates - so that 
the term on the right of (119) will drop out - we now come to the less trivial category 
of applications for which the term "virial theorem" is most commonly employed [9], 
namely cases for which the weighting factor X has homogeneous quadratic dependence 
on the Cartesian space coordinates X\ The simplest such possibility is the isotropic 
case obtained by taking X proportional to the square of the radial distance r from the 
center as given (in a Cartesian system with origin at the center of mass) by 

= ^ijX'X^ (128) 

while another obviously important special case is that for which X is taken to be 
proportional to the square of the distance zu from let us say the axis, as defined by 

= (7,,- - 61S])X^X^ = {X^y + (X^)2 . (129) 

The isotropic case obtained simply by identifying X with is that for which A4 is 
just the ordinary scalar quadrupole moment /. Thus by setting 

X^r"^ =^ M^Q, WiWjX ^ , (130) 

one sees from the generic virial relation (119) that the scalar quadrupole moment 

Q = J pr^d^X. (131) 

will satisfy a time evolution equation of the form 

^^3 = 2/r;,d3x. (132) 

It evidently follows that for a stationary (i.e. time independent) configuration the 
integral on the right of this equation must vanish. 

The axial case obtained by identifying X with w"^ is that for which A4 is just the 
moment of inertia moment about the X^ axis. Thus by setting 

X^w^ =^ M^I,, ViVjX ^ 2{-fij - S\5]) , (133) 

one sees from the generic virial relation (119) that this moment of inertia 

4 = y pw^d^X, (134) 

will satisfy a time evolution equation of the form 

^^3 = 2/(rj/ + rf)d^x. (135) 
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It evidently follows again that for a stationary (i.e. time independent) configuration 
the integral on the right of this equation must vanish. 

The preceeding special isotropic and axial cases can be considered as combinations 
of the separate components of the corresponding tensorial theorem - as obtained by 
substituting X'^X^ in place of A" in (117) and (119) - according to which the first and 
second time derivatives of the generic quadrupole moment component 

g^i = j pX'X^d^X, (136) 

will satisfy equations of the form 

^Q'^ = 2 I n^'x^^ d?X 
at J 

and 

^,Q'' = 2 / d^X . (138) 

As well as the stationary case, for which the integrals on the right of (137) and (138) 
must vanish, the time dependent case is also of particular interest because it is the third 
time derivative of the trace free part of the quadrupole moment tensor that provides 
the source of gravitational radiation in the Newtonian (weak field, low velocity) limit 
of general relativity. 

The foregoing conclusions, and in particular the scalar quadrupole evolution equa- 
tion (132), can be construed as a concise general statement of what is commonly re- 
ferred to - in more detailed applications to particular well known cases - as the "virial 
theorem". The most widely familiar version [9] applies just to the case of a single con- 
stituent perfect fiuid, but in a recent study of ellipsoidal configurations Sedrakian and 
Wasserman [10] have provided a version applying to cases in which there is a pair of 
independently moving fluid constituents. Our present - formally very simple - result 
(132) actually goes further, not just by allowing for the possibility of more than two 
constituents, but also, less trivially, by including allowance for the effect of entrainment 
whereby the strong coupling between the constituents modifies the momenta of the con- 
stituent particles whose effective masses may deviate substantially from their ordinary 
rest masses. 

To relate our concise new general purpose "virial theorem" (132) to the more spe- 
cialised results that are already well known, it is instructive to consider the distinct 
contributions that are involved. To start with it can be seen from (38) that the total 
energy (106) will be expressible in the form 

^tot = ^.at + ^.H , (139) 

where the purely material contribution is defined by 

^.at= lu^..d'X, (140) 
22 



(137) 



and the total gravitational binding energy contribution is defined by 

= E^, + E^^, = j U^^, d'X , (141) 



grt 



in which 



grf 



and 



E„ = ju„'^'X = j ■PI'A'^- (143) 

For a confined source, it can be seen from (40) using Green's theorem that the to- 
tal gravitational contribution (141) will be related to the separate gravitational field 
contribution (142) and gravitational potential contribution (143) by 

^g. = -E^, = l^pot ■ (144) 
In a similar manner, using the observation that 

T^^i - -U,. ' (145) 

it can be seen that the integral on the right hand side of (119) can be split as a sum of 
a purely material contribution and a gravitational contribution, in which (again using 
Green's theorem) the latter works out to be the same - for a confined source - as the 
corresponding contribution to the energy, i.e. one obtains 

JT:^d'X = JT:^.^^"X + E,... (146) 
The material energy contribution can evidently be further decomposed as 

^.at = ^Mn + E,^, , (147) 

where the separate kinetic and internal contributions are expressible in terms of the 
densities given by (18) as 

£^.n = / ^Mn , E,^^ = J U,^, d^X . (148) 

of which the latter will be expressible in terms of the rest frame chemical potentials 
introduced in (19) as 

E,^, = E / n,x" d?X- jm d?X . (149) 
The corresponding decomposition for the material contribution in (145) has the form 

/ ^.L. = / TV d'X + / TV d'X , (150) 



23 



in which the kinetic contribution can be seen from (15) and (18) to be given by 

jTJ^d'X = 2E^^, (151) 
while the internal contribution can be seen from (16) to be given by 

J Tj, d=^X = 3 y * d^X + E / Kxl . (152) 

Putting the separate pieces together, it can be seen that the integral in the for- 
mulation (132) of the generalised "scalar virial theorem" - and that in any stationary 
equilibrium configuration must therefore vanish, i.e. 

jT:,.d'X = 0, (153) 

- will be expressible by 

/ TJ, d'X = E^, + 2E^^ + ^j^d'X + Y. j Kxl d'X . (154) 

In the combination on the right hand side - which must vanish for a stationary equilib- 
rium configuration - it is to be remarked that the first two terms (the total gravitational 
energy plus twice the ordinary kinetic energy) are of the kind that is traditionally fa- 
miliar. As compared with the multiconstituent version given recently by Sedrakian 
and Wasserman [10] for the idealised case in which mutual entrainment between the 
constituents was ignored, the difference here consists firstly of the use of the undecom- 
posable pressure function ^ in place of a sum (of the form E-Px) ordinary pressure 
contributions from the distinct constituents, and secondly of the inclusion of the final 
term involving the previously ignored entrainment momenta themselves. 

Before concluding, we wish to point out that the homogeneously quadratic category 
that has been discussed in detail immediately above is not the only useful category 
of non-trivial applications of our generic virial theorem (119). Another category worth 
consideration is that of homogeneously first (as opposed to second) order mass moments. 
Strictly linear dependence on the Cartesian space coordinates leads merely to the trivial 
dipole moment case discussed above, but a no less natural, and much more interesting, 
alternative possibility is that of the isotropic homogeneously first order (but non-linear) 
case obtained by setting X — while the next most obviously interesting possibility 
is that of the cylindrical homogeneously first order (but non linear) case obtained by 
setting X = w. 

For the first of these, namely the isotropic first order case, we have 

X =^ ViVjX = r~-^(7jj - VriVrj) , (155) 

where is the radial unit vector as defined by 

i/; = r-^X' . (156) 
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Since, in terms of standard spherical coordinates r, 6, ip, the volume element will be 
given by dX^ dX^ dX^ = sin 9 dr d^^ d(/9, the isotropic homogeneously linear virial 
theorem obtained by substituting (155) in (119) will be expressible as 

— j pr^smOdrddd^p^ j (T ^ - Tji/^iZ/^j)rsin^drd^d(/?. (157) 

An even simpler relation is obtainable for the cylindrical first order case, for which 
we have 

X = w =^ ViVjX = w~^u^piU,^j , (158) 

where is the unit vector in the direction of the relevant axial Killing vector (60), 
namely 

u;^w-^k^, (159) 

so that its components will be given by = —X'^/zu, = /w, = 0. This leads 
to what is interpretable as the Newtonian limit [11] of a result originally derived in a 
relativistic context by Bonazzola [12] . In terms of standard cylindrical coordinates, as 
specified by X^ = lucosip, X^ = w sirup, X'^ = z, the volume element will be given by 
dX^ dX^ dX^ = w dw dz dp. It can thus be seen that the cylindrical homogeneously 
linear virial theorem obtained by substituting (158) in the generic relation (119) will 
take the form ^ 

J pw"^ dw dzdp = j T^^^u^piU^j dw dz d(p . (160) 

More particularly, in the stationary axisymmetric case, for which the left hand side 
vanishes and for which the ip integration is trivial, (160) will reduce to the 2-dimensional 
integral relation 

'rju^.u^jdwdz^O, (161) 



I' 



that was originally derived by Gourgoulhon and Bonazzola [11], in which it can be seen 
from (29) that the gravitational contribution to the integrand can be separated out in 
a decomposition that takes the form 

Ty^^^^J = T^Jl^;l^^-U^,, (162) 

as a consequence of the axisymmetry. 

We conclude by emphasising that, like the ordinary virial theorem (153), the Bonaz- 
zola theorem (161) is a particular example of the condition 



I' 



TjyiWjX d^X ^ , (163) 

that can be seen from (119) to be a generic requirement for any stationary equilibrium 
configuration for all admissible choices of the weighting factor X. It is to be recalled that 
T^*-^ is the total of all material and gravitational stress contributions, and that, for the 
purpose of this generic theorem, the relevant admissibility criterion is that X should be 
any given function of the Cartesian space coordinates whose gradient components VjA" 
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increase more slowly than a quadratic function of these coordinates at large distance. 
(This condition would fail for a substitution of the form X r^, but it is evidently 
satisfied by the substitution X that leads to the ordinary virial theorem, and 

a fortiori by the substitution X ^ zu that leads to the Bonazzola virial theorem.) 
Whereas the individual terms in the expansion (152) have forms that depend on the 
specific nature of the multiconstituent fluid models developped above, on the other 
hand the generic relation (163) and its dynamical generalisation (119) depend only on 
the generic form of the energy momentum and mass conservation laws (110), (111), 
and (113), which should hold for any (complete) Newtonian continuum model. These 
generally applicable laws are all that is needed for the derivation of the local virial 
equation (114) that underlies the global relation (119). The generic virial equilibrium 
relation (163) is just the global consequence of the local equilibrium condition 

V,V,7;:f = 0, (164) 

that obviously holds as the stationary specialisation of the local virial equation, namely 

the dynamical relation (114), whose importance - as an easily memorable law that 
must be satisfied by any Newtonian continuum configuration - has not been adequately 
recognised hitherto. 
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Appendix: Noether identities in Newtonian theory 

Let us consider the generic case of a Newtonian model characterised by a total action 
density that is a sum 

a 

of contributions labelled by some index a. This includes the kind of model dealt with 
in the preceeding work which can be described by taking a to range over the four 
values specified by the labels "kin", "int", "pot" and "gra". In what follows we shall 
use the symbol = to denote equivalence modulo a space time divergence. The content 
of the action principle is thereby expressible as the postulate that when the relevant 
field equations are satisfied, any "admissible" infinitesimal variation of the relevant 
dynamical field variables will give rise to a corresponding total variation ^A^^^ that is 
equivalent to zero in this sense, i.e. ^A^^^ = 0, so that its space time integral will vanish 
by Green's theorem for any variation with compact support. 

The purpose of the present section is to consider the effects of more general - not 
necessarily "admissible" - variations of the relevant dynamical fields as well as variations 
of the various spacetime background fields on which the complete specification of the 
action depends. For such more general variations, the total action density variation will 
satisfy a relation of the form 

<^Atot = E^*A„ + ^5«A„, (166) 
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in which denotes the contribution from the variations of the background fields and 
5" denotes the contribution from any non-admissible parts of the variations of the dy- 
namical fields. 

In the technically simpler case of a relativistic model the relevant spacetime back- 
ground variation would be fully determined by the variation of the space time metric 
Qij^i,, or equivalently of its contravariant inverse g^^". However, in the Newtonian case 
under consideration here, it is not sufficient to know the variation of the corresponding 
contravariant space metric 7^^. It will be necessary as well to know the associated 
variation of the corresponding preferred time gradient t^. To deal with more elaborate 
models, the variation of the linear connection might also be involved, but this will not 
be necessary for models of the simple kind considered here. However to deal with the 
Galilean gauge dependent contribution A^.^ it will also be necessary to take into account 
the variation of the chosen ether frame vector e'^. Thus the relevant background field 
variations will be given by expressions of the form 

in which the final term will drop out for gauge independent contributions such as A^^^ , 
Ap^j, and A^^j. In the models under consideration, the relevant dynamical fields are the 
gravitational potential, 0, for which arbitrary variations are admissible, and the current 
vectors whose admissible variations are restricted, so that a generic variation Sn^ 
may include a non - admissible part S'^n^ that provides a contribution of the form 

^'^^-^/<- « 

As in the (for this purpose simpler) relativistic standard procedure [13] for 

the derivation of useful Noether type identities is to consider variations of the trivial 
kind generated by an arbitrary displacement vector field, .^^ say. This means that the 
variation of each (background or dynamical) field variable will be given by the negative 
of its Lie derivative. The relevant formulae are thus given by 

-6Aa = i£Aa^e'^p^a, (169) 

-S-f"" = i£ = ^"^pl""" - 27" ^'^ VpC"^ , (170) 

-5t^ = i£, = ^"Vpt^ + ipV^e^' , (171) 

-Se^' = i£. = eV.e'^ - e^'V^^ , (172) 

and 

-5< = i£ < ^ e'^.n^ - <Vpe ■ (173) 

In view of the uniformity properties (1-16) and (1-18) of the unperturbed background 
fields, we shall be left with 

S-fi"' = 2Y^''VpC^ , (174) 
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St^^-t,V^C', (175) 

Se" = e^VpC" , (176) 

and it can be seen from the expression (1-86) for the form of an admissible current 
variation that the only remaining (i.e. inadmissible) contribution from (173) will be 
given simply by 

5«< = <Vpe- (177) 

Since (using the notation introduced in (166) for equivalence modulo a divergence) 
the variation (169) will evidently satisfy 

5Aa ^ AaV^e (178) 

it can be seen that the terms in (166) can be regrouped on the left to give a relation of 
the simple form 

2:„r.v^r = o, (179) 

in which the quantity that will be interpretable as the relevant total stress momentum 
energy density tensor can be read out in the form 

TJ. = j:Ta'., (180) 

a 

as the sum of contributions given by the formula 

^ ^ (5^ - 2 + ^ t -^^e^ (181) 

in which the relevant generalised pressure contribution ^'^ is given in terms of the 
corresponding momentum contributions, 



X 



by 

It is to be remarked that, due to the restrictions 

= -7^" 5U , U Se" = -e" SU , (184) 

resulting from (1-3) and (1-16), there is some gauge ambiguity in the specification of the 
partial derivative coefficients introduced in (167) but it can easily be checked that the 
stress momentum energy contributions T^^, specified by the combination (181) will be 
physically well defined in the sense of being unaffected by the choice adopted in (166). 

By a further equivalence transformation (modulo a divergence) it can be seen that 
(179) can be converted to the form 

rv,3:„r.-o. (i85) 
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Since this must hold for an arbitrary vector field and hence in particular for a 
displacement field with compact support in any small space-time neighbourhood, it can 
be seen by integrating over such a neighbourhood (so that the divergence ambiguity in 
the equivalence relation cancels out by Green's theorem, as in the usual derivation of 
the dynamical equations from the action principle) that the coefficient of in (185) 
must vanish, i.e. we obtain a total energy momentum conservation law of the form 

VMrr. = 0. (186) 

It remains to show that the conserved total stress momentum energy density ten- 
sor T^l'i, obtained by the foregoing procedure is actually the same as the tensor T^^^i, 
introduced in (30) 

The relation (186) is not an identity in the strict sense, since its derivation depends 
on the dynamical field equations obtained from the action principle. To obtain a purely 
mathematical identity it is necessary to take account of the unrestrained variations of 
all relevant field variables. In the application with which we are concerned here this 
includes not just the background field variations involved in (167) and the variations 
of the currents but also the variations of the gravitational potential (j) and its gra- 
dient, so that the complete generic expression for the variation of an action density 
contribution will have the form 



(9A 8\ 

5Aa = S^K + E K + ^50 + ^ 5(0,,) . (187) 

For variations generated, as before, by an arbitrary displacement vector field we can 
again use the background field variation formulae (174), (175) and (176) together with 
the corresponding dynamical field variation formulae (173) and 

50 = -e£0 = -e0,,, (188) 

we thus obtain a relation of the form 

which will be satisfied as an identity in the sense of being independent of the dynamical 
field equations. Since it is possible to choose both and its derivative V,^*^ indepen- 
dently at any given point, it follows that the corresponding coefficients must vanish. 
Thus from the left hand side of (189) we obtain the obvious identity 

(9A (9A 

V.A„ ^TK, V.n^ + — ^ + — -\ V,V.0 , (190) 
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while from the righthand side we obtain an identity of the less trivial form 

- E n^\K - 7- + ^ - 1^ e'^ - . (191) 



A noteworthy special case is that of the internal contribution A.^^., which depends 
neither on the gravitational potential (unlike Ap^J nor on the ether frame vector 
(unhke Aj^.^), so that when it is substituted for A^ in (191) the first and the last term 
will both drop out. Contraction with and 7'^'^ can then be used to eliminate the other 
terms at the end, leaving a result that can be recognised as the simple identity (1-143) 
that was quoted above, while the following identity (1-144) can similarly be derived by 
performing an antisymmetrisation after contracting just with 7^^^. 

Quite generally, the Noether identity (191) can be used to replace the formula (181) 
by the more explicit and manifestly gauge independent expression 

dA 

r„t = ^aS't + E vr-. < - . (192) 

With the aid of (190) it can be directly verified that each such contribution will satisfy 
a divergence identity of the form 

V^T„t = E/„^. + ^V.0, (193) 

in which the Eulerian derivative has a definition of the usual form 

while the force density contributions are given by an expression whose form 

f^, = 2<V[.7r,%] + TTj^V.^ . (195) 

is analogous to that for the combined force densities (1-159). 

By summing over these separate contributions we obtain the corresponding formula 

= ^.Jt + E ,n-^K - 1^ , (196) 

for the total stress momentum energy density tensor, where the corresponding total 
pressure function and momenta are given by 

a a 

while there will be an analogous expression 

/J. = E/a^' (198) 
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for the total non-gravitational force density exerted on each constituent. 

In the application considered in the preceeding work the purely gravitational action 
density A^^j given by (25) provides no contribution to the momenta, i.e. we shall 
have TT = 0, and hence ^ , = A , so we can make the identifications tt — tt^, 

grf ' grf grf tot ' 

f ^„ — ff,, and ^f. . = + A , in which it is to be recalled that we use a blank label 

tot ^ J 1/ ^ tot grf 

to indicate the sum over all values of a except the gravitational contribution "gra", 
i.e. for the sum over the values "kin" , "int" , and "pot" , which is the only part that is 
relevant if (as in the Cowling approximation) one is concerned just with evolution in 
a fixed gravitational background, but not with the effects of self gravity. It can thus 
be verified that the definitions introduced in the systematic mathematical procedure 
developed in this appendix arc entirely consistent with those introduced ad hoc in the 
main part of the text. In particular, it can easily be seen from (196) that the general 
purpose prescription provided by (181) and (182) for T^^^ leads to a result that is in 
exact agreement with what is given by the formula (30) that was obtained from more 
specific physical considerations in the main part of the text, while similarly by summing 
over the index a in the identity (195) one recovers a result that can be seen to agree 
with the previously quoted divergence formula (3). 
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